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ABSTRACT 

We study close encounters involving massive main sequence stars and the evolution 
of the exotic products of these encounters as common-envelope systems or possible 
hypernova progenitors. We show that parabolic encounters between low- and high- 
mass stars and between two high-mass stars with small periastrons result in mergers 
on timescales of a few tens of stellar freefall times (a few tens of hours). We show 
that such mergers of unevolved low-mass stars with evolved high-mass stars result in 
little mass loss (~ 0.01 M Q ) and can deliver sufficient fresh hydrogen to the core of 
the collision product to allow the collision product to burn for several million years. 
We find that grazing encounters enter a common-envelope phase which may expel 
the envelope of the merger product. The deposition of energy in the envelopes of our 
merger products causes them to swell by factors of ~ 100. If these remnants exist in 
very densely-populated environments (n > 10 7 pc -3 ), they will suffer further collisions 
which may drive off their envelopes, leaving behind hard binaries. We show that the 
products of collisions have cores rotating sufficiently rapidly to make them candidate 
hypernova/gamma-ray burst progenitors and that ~ 0.1% of massive stars may suffer 
collisions, sufficient for such events to contribute significantly to the observed rates of 
hypernovae and gamma-ray bursts. 
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1 INTRODUCTION 

Stellar collisions are an important physical process. Al- 
though stars are very small in comparison to the typical 
sizes of stellar clusters, long stellar lifetimes, high number 
densities, gravitational focusing, or a combination of these 
factors can permit a significant fraction of the populations 
of some stellar systems to undergo collisions. 

Globular clusters, with number densities of ~ 10 4 
pc -3 , velocity dispersions of ~ 10 km s -1 and long — lived 
low-mass stellar popul ations, have long been known to 
host stellar collisions. iHills fc Davl <ll976l) estimate that 
40% of the stars in some globular clusters have suffered 
collisions. The remnants are thought to appear on cluster 
colour-magnitude diagrams as part of the cluster's popu- 
lation of blue stragglers, although it does not appear that 
all b lue stragglers are products of collisions jDavies et alJ 
120041) . The collisional formation of these object s has been 
studie d hydrodynamical l v bv| Benz fc Hills) |l987h.|Lai e t alJ 
Jl99ot.lLombardi et ail <2002lh ISandauist et alflll997l) and 
ISills et all tOOti) . 
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Work bv lBonnell et al.l (Il998l) showed that star clusters 
initially experience a phase of contraction which can drive 
the stellar density in the core to extreme values (~ 10 8 
pc -3 ). Under these conditions, encounter timescales can be 
sufficiently short to allow short-lived high-mass stars to ex- 
perience collisions. Older stellar systems with high (but not 
extreme) number densities, such as Orion, will also rarely 
host collisions involving massive stars. The appearance and 
properties of merger remnants depend cruci ally upon the 
degree of mixing induc ed by th e collision. |Benz fc Hill " 
<1987jh iLombardi et al] j2002lh ISandauist et alJ dl99 ' 
and [silk^tHji] I2002T) all studied encounters between 
low-mass ma i n seq uence stars and, with the exception of 
iBenz fc Hilll (1987), found that mixing was minimal. It is 
important to determine if this is also the case in encounters 
involving high-mass stars, since the merger may deliver 
fresh hydrogen to the core of the massive star, potentially 
reju venating it. 

iPorteeies Zwart fc McMillan! (I2002T) proposed that the 
extreme densities thought to exist in some very young 
stellar systems may trigger a runaway merger process, 
leading to the formation of an intermediate-mass black 
hole. These objects have been posited as an explanation 
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for U ltra-Luminous X-ray sources (e.g. iTaniguchi et alJ 
(2000)) and as the building blocks for the supermas- 
sive black holes res iding in the cores of many galaxies. 
iFreitae et all J2005al) conducted Monte Carlo simulations of 
star clusters experiencing core collapse and hence achieving 
sufficient stellar densities for runaway collisions to occur. 
To include the effects of ste l lar co llisions in stellar cluster 
simulations, IFreitae fc Bend i2005l) performed a very large 
suite of smoothed-particle hydrodynamics (SPH) collision 
simulations covering a large range of stellar masses (0.1 — 75 
Mq) and collision velocities. 

In this work, w e study a sma l ler pa rameter space 
than that covered bv IFreitae fc Beirzl i2005h . but our SPH 
calculations are at higher resolution (typically involving 
~ 7 x 10 4 particles) and include the highest-resolution 
simulations of massive stellar encounters ever performed, 
using ~ 6 x 10 5 particles. The use of two very different 
resolutions allowed us to demonstrate that our conclusions, 
in particular the rapid mergers we observe in encounters 
with small periastrons, are independent of numerical effects. 
Our superior resolution also allowed us to examine aspects 
of collisions and encounters between h igh-mass main 
sequence stars that IFreitae fc Bend <l2005f) did not study, 
such as the rotation of the collision products. We also 
performed simulations in which the enco unter periastron 
was la rger than the sum of the stellar radii - iFreitag fc Benzl 
(2005) explicitly neglected such encounters. Our work is 
particularly relevant to determining the properties and 
evolutio n of collision products and to t he ru naway mergers 
model of IPortegies Zwart fc McMillanl (120021) . 

In Section 2, we discuss our numerical techniques and 
the calculations we have performed. In Section 3 we give 
the results of our calculations. In Section 4 we discuss the 
possible outcomes of encounters with the collision products 
described in Section 3. We present our conclusions in 
Section 5. 



2 NUMERICAL TECHNIQUES 

The SPH code used in this work is fully described in 
i Batd J1995D. The cod e uses the cubic spline kernel of 
iMonaghan fc Lattanziol il985l) and enforces a neighbour 
number of ~ 50 fo r all particles. We used t he a rtificial viscos- 
ity prescription of lOingold fc Monagharj (|l983f ) , with a = 1 
and f3 = 2. The majority of our SPH calculations were com- 
posed of ~ 70, 000 particles and were run on Linux PCs. We 
also performed two very high-resolution calculations com- 
posed of « 623, 000 particles on the SGI Origin 3800 at the 
United Kingdom Astrophysical Fluids Facility (UKAFF). 

We generated spherically-symmetric equilibrium stel- 
lar models us i ng th e ATON stellar evolution code (e.g. 
IVentura et alJ l)l998h '). We extracted from these models ra- 
dial profiles of density, internal energy and helium frac- 
tion and constructed three-dimensional (although initially 
spherically-symmetric) SPH models of the stars. 

We utilised two solar-metallicity stellar models - a 1 
Mq star approximately halfway through its main sequence 
life, and a 9 Mq star near the end of its main-sequence 
phase. We constructed our SPH stellar models by building a 
uniform close-packed grid of equal-mass SPH particles and 



iteratively adjusting the mass of each particle until the dis- 
crepancy between the local SPH density at a given radius 
and the corresponding density derived from the ATON code 
was less than 0.1%. This procedure was able to accurately 
model the 1 Mq star using 11,000 SPH particles and 
the 9 M Q star using > 611,000 particles for the UKAFF 
calculations without generating any features pathological to 
the simulations. However, the density contrast between the 
core and the envelope of the 9 Mq star was so great that 
it could not be represented by a uniform grid of ^ 60, 000 
particles in our PC runs. We observed that the envelope of 
the object was poorly resolved and that the iteration proce- 
dure had generated a population of very low-mass particles 
in the core, a situation which can lead to unphysical re- 
sults. This was a consequence of the stellar core, in which 
the density changes sharply over a short lengthscale, be- 
ing only a few smoothing lengths across. Our solution was 
to build a uniform-density grid, containing a large enough 
space-density of SPH particles to resolve the core, out to a 
radius r gr id < r max , where r max is the 'real' radius of the 
star (see below), and to deform the grid outside the core 
radius r core (taken to be 1 Rq) stretching it to a radius of 
r m ax, leaving a uniform particle distribution inside the core 
while redistributing particles in the envelope to allow better 
resolution near the surface of the star. The deformation was 
radial, achieved by multiplying the x, y and z components 
of each particle's initial position vector r by a factor / given 
by 

{1 for |r| ^ r core 

i+^ferfe- 1 ) fOT r grld > \r\ > r core 

Cumulative mass and density profiles for the two SPH 
model stars are shown in Figure Q and compared to the 
one-dimensional models from the ATON code. Figure^also 
shows the helium fraction profile of the stars. The fits of 
the three-dimensional SPH models to the one-dimensional 
ATON model are clearly very good. Note that none of the 
9 Mq models has a mass of exactly 9 Mq. The mass of the 
ATON model star is 8.85 Mq, that of the high-resolution 
SPH model is 8.75 Mq and that of the low-resolution 
model is 8.56 Mq. The SPH models have lower masses than 
the ATON model because we chose to neglect material of 
extremely low density beyond a radius of ~ 5 Rq in the 
ATON model. The radius at which this cutoff occurred 
resulted in the UKAFF and PC models having slightly 
different masses. The radius of an SPH stellar model is the 
radius out to which the particle distribution extends plus 
the average smoothing length of the particles at that radius. 
Since the particles in the high-resolution model had smaller 
smoothing lengths, more particles could be fitted inside 
the radial cutoff, leading to a higher mass. We allowed all 
the SPH models to relax before beginning any simulations. 
Both massive SPH model stars expanded on relaxation, 
the low-resolution model slightly more so because of its 
lower mass, leaving the massive stellar models with radii of 
~ 6 Rq. The convection zone delimiting the core of the 9 
Mq star is visible in Figure as the discontinuities in the 
cumulative mass and density profiles at a radius of w 1 Rq . 
The cumulative mass plots of the 9 Mq star imply that the 
core mass of this star is « 2 Mq and the helium fraction 
plot shows that the 9 Mq star has almost exhausted its 
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Figure 1. Cumulative mass, density and composition profiles for the two main-sequence stellar models used in our calculations (dot- 
dashed lines: relaxed high-resolution UKAFF models, dashed lines: relaxed lower-resolution PC models, solid-lines: ATON models). 



Run no. 


Mi (M ) 


M 2 (M ) 


Pcriastron (R©) 


r P /(Ri - 


- R2) No. particles 


Remarks 


1 


8.56 


1.00 


1.0 


0.13 


68, 000 




2 


8.56 


1.00 


2.0 


0.26 


68, 000 




3 


8.56 


1.00 


3.0 


0.39 


68, 000 




4 


8.56 


1.00 


4.0 


0.51 


68, 000 




5 


8.56 


1.00 


5.0 


0.64 


68, 000 




6 


8.56 


1.00 


7.0 


0.81 


68, 000 




7 


8.56 


1.00 


10.0 


1.30 


68, 000 




8 


8.56 


1.00 


15.0 


1.95 


68, 000 




9 


8.56 


1.00 


2.0 


0.26 


57, 000 


as 2 using a point mass 1 Mg star with a 
gravitational smoothing length of 0.5 Rq 


10 


8.56 


1.00 


2.0 


0.26 


57, 000 


as 2 using a point mass 1 Mg star with a 
gravitational smoothing length of 0.05 Rq 


11 


8.76 


1.00 


1.0 


0.13 


623, 000 


as 2 but higher resolution 


12 


8.76 


1.00 


2.0 


0.26 


623, 000 


as 2 but higher resolution 


13 


8.56 


8.56 


1.0 


0.07 


114, 000 




14 


8.56 


8.56 


2.0 


0.14 


114, 000 




15 


8.56 


8.56 


3.0 


0.21 


114, 000 




16 


8.56 


8.56 


4.0 


0.28 


114, 000 




17 


8.56 


8.56 


5.0 


0.35 


114, 000 




18 


8.56 


8.56 


10.0 


0.70 


114, 000 




19 


8.56 


8.56 


15.0 


1.05 


114, 000 




20 


8.56 


8.56 


18.0 


1.26 


114, 000 





Table 1. Runs performed. Each run is defined by the masses of the stars, the periastron at which the encounter occurs and the total 
number of SPH particles used in the simulation. 
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core hydrogen. The SPH model 1 Mq stars also expanded, 
but the cumulative mass and density plots show that the 
agreement between the SPH models and the ATON models 
is still very good throughout almost all of the star's volume. 
The core of the 1 Mq star, delimited by the volume of the 
star within which hydrogen burning is in progress, has a 
radius of as 0.2 Rq and a mass of ~ 0.3 Mq. Both stars, 
particularly the 9 Mq star, are strongly centrally condensed 
and the 1 Mq star is considerably denser throughout almost 
all of its volume than even the core of the 9 Mq star. 

Once relaxed, the model stars were placed on parabolic 
orbits with a variety of periastrons and sufficiently far apart 
that there was negligible tidal interaction at the beginning 
of each simulation. 

Encounters occurring at periastrons r p < Ri + R2, 
where Ri and R2 are the radii of the stars, necessarily result 
in a physical collision on the first periastron approach. If 
the two stars do not suffer such an immediate collision, 
they raise tides in each other and some of their orbital 
kinetic energy is converted into thermal and kinetic energy 
in the stars' envelopes. If the transfer of orbital energy to 
tidal energy is large enough to leave the total orbital energy 
negative, the stars become bound to each other. The energy 
required to bind the stars is dependent on their relative 
velocity at infinity Voo- 

All encounters described in this work are parabolic, 
so that Voo — 0. However, from the orbital energy at the 
end of each of our simulations, we calculate the value of 
Voo that would have rendered the system just bound, which 
we term Avoo- We can thus determine which encounters 
would result in bound systems in any cluster environment, 
provided that Voo is small compared with the relevant 
stellar escape velocities, while retaining the computational 
simplicity of performing parabolic calculations. This paper 
is concerned with stellar encounters in clusters with low 
velocity dispersions, ~ 10 km s _1 . This sets the 'energy 
scale' of our simulations. 

The physical parameters defining each of our runs are 
the stellar masses Mi, M2 and the orbital periastron r p . 
We list all runs performed in Tabled 



3 RESULTS 

In Figure we plot the evolution of the thermal, kinetic, 
gravitational potential and total energies during the high- 
resolution UKAFF calculation involving a 9 Mq star and a 
1 Mq star with a periastron of 2 Rq. In Figure [3J we plot 
the separation of the stellar cores (taken to be equal to the 
separation of the stars' centres of mass) against time for 
the same simulation and compare it with the corresponding 
PC run. We take this as our 'fiducial run', against which 
we compare all other runs. 

There are 4 distinct peaks visible in the energy and 
core separation plots corresponding to 4 periastron passages 
occurring before the system formed a single merged object. 

In Figure [I] we follow the mass-loss in this calcula- 
tion. We assume that any SPH particles with positive total 
energy and moving away from the centre of mass will be 
lost from the system. Figure 0] shows that mass ejection 
in this simulation proceeds in several distinct phases, but 
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Figure 2. Evolution of system energy components during the 
high-resolution encounter between a 9 Mq star and a 1 Mq star 
with a periastron of 2 Rq. The corresponding plots for the the 
low— resolution calculation are not shown because they are difficult 
to distinguish. 



that the total quantity of mass lost by the time the cores 
merge is not large - a few tenths of a solar mass or a few 
percent of the combined stellar mass. Figures |3] and 2] show 
that the mass-loss rate increases sharply as the intruder 
star approaches its first periastron and decreases after the 
first periastron approach. The rate falls almost to zero at a 
time of ~ 10 hr because, as Figure [3] shows, the 1 Mq star 
leaves the 9 Mq star's envelope. At a time of ~ 12 hr, the 
intruder falls back into the larger star and mass ejection 
begins again. The 1 Mq star has now lost enough kinetic 
energy that it is unable to leave the envelope of the 9 Mq 
star again and spirals rapidly towards the larger star's core. 
The cores merge at « 18 hr, when the core separation 
becomes effectively zero. Mass loss ceases rapidly after this 
point. 

Quantitative differences between the UKAFF and PC 
runs were very small. The slightly greater mass of the 9 
Mq star in the UKAFF runs led to slightly tighter orbits, 
greater dissipation of kinetic energy (by sa 1%) and to 
mass loss smaller by a few hundredths of a solar mass. 
We conclude from the concordance of the results that the 
high-mass stars studied in this paper can be accurately 



modelled using ~ 6 x 10 particles. 



3.1 Capture radius 

To calculate the distance within which the stars must 
approach to become bound to each other for a given 
value of Voo, we determined the orbital parameters of each 
simulation after the first periastron passage. We used an 
iterative procedure to decide which (if either) star each SPH 
particle was bound to and calculated the centre-of-mass 
positions and velocities of both objects. From these, we 
determined the total orbital energy and hence Avoo for 
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5 10 15 20 25 30 
Time (hours) 

Figure 3. Separation of stellar cores during the high-resolution 
(solid line) and low-resolution (dashed line) encounter between a 
9 Mq star and a 1 Mq star with a periastron of 2 Rq. 



each initial encounter. In Figures |S] and HJ we plot Av^ 
as a function of periastron for encounters between 9 Mq 
stars and 1 Mq stars and encounters between two 9 Mq 
stars respectively. Note that the last point on each plot 
should be treated with some caution, since the error in 
the value (determined by the accuracy of the SPH code's 
energy conservation) is comparable to the measurement. 
The error in the code's energy conservation is fractional 
- in both cases we have multiplied the fractional error by 
the total energy of each run (which is different for the two 
types of encounter considered) and converted the result 
into a velocity. We are interested in stellar systems with 
low velocity dispersions such as the Trapezium cluster, so 
we define the capture radius, r c , to be the periastron at 
which Avoo = 10 km s _1 . Note that we have here implicitly 
assumed that Avoo is independent of v x . We repeated the 
encounter between two 9 Mq stars with a periastron of 
15J?q on three hyperbolic orbits with Voo — 10, 20, and 50 
km s _1 and found that this assumption is justified within 
this range of velocities, which covers the stellar systems in 
which we are interested. 

Figures and |S| show that the capture radius for 
encounters between 9 Mq stars and 1 Mq stars is w 11 
Rq and that for the encounters between two 9 Mq stars 
is ~ 18 Rq. The radius of the 9 Mq star is m 6 Rq, 
implying that encounters between 9 Mq stars and 1 Mq 
stars with periastrons between 7 and 11 Rq and encounters 
between two 9 Mq stars with periastrons between 12 and 
18 Rq lead initially to a bound system of two stars. The 
borderline between what is considered a 'captured system' 
and an 'immediate merger' is somewhat arbitrary. For 
example, encounters between 9 Mq stars and 1 Mq stars 
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Figure 4. Mass loss during the high-resolution (solid line) and 
low— resolution (dashed line) encounters between a 9 Mq star and 
a 1 Mq star with a periastron of 2 Rq. 
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Figure 5. Plot of Ad^o after first periastron passage against ini- 
tial periastron distance for encounters between 9 Mq stars and 1 
Mq stars. The dotted line is a fit to the measured values and the 
solid line represents Avcc = 10 km s — 1 . The error on each point 
is ±8 km s —1 . 

with initial periastrons of 4 and 5 Rq, which are strictly 
grazing collisions, produced captured systems in which the 
intruder stars were able to exit the envelopes of the 9 Mq 
stars on very eccentric elliptical orbits. Such systems cannot 
be studied further by means of SPH simulations because 
the CPU time required is prohibitively large. However, we 
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Figure 6. Plot of Auoo after first periastron passage against ini- 
tial periastron distance for encounters between two 9 Mq stars. 
The dotted line is a fit to the measured values and the solid line 
represents Avoo = 10 km s — 1 . The error on each point is ±17 km 



can study their evolution indirectly, as described in the 
next section. 



3.2 Evolution of grazing encounters 

In encounters where the intruder star passes within the cap- 
ture radius but suffers only a grazing collision or misses the 
other star completely, a captured system consisting of two 
stars in highly eccentric orbits results. Although we cannot 
study the evolution of these systems hydrodynamically, we 
can attempt to infer their evolution from the orbits in which 
they are left by their first encounter. 

We first assume that, since the periastron changes very 
little in tidal encounters, the quantities of tidal energy AE P 
and angular momentum AJ P extracted from the orbit on 
the first periastron passage are also extracted at each sub- 
sequent periastron passage, and that the tidal energy taken 
up by the stars themselves is efficiently radiated away. For 
the nth periastron passage, we therefore have, for the energy 
E n and orbital angular momentum J n 



En — Einitial TlAEp 



-nAE v 



(2) 



where Einitial is the initial orbital energy and is zero for 
parabolic encounters, and 



Jn — Jit 



nAJp, 



(3) 



where Jinitial is the initial orbital angular momentum. 

We then calculate the eccentricity e„, semi-major axis 
On, apastron r n and period T n for each orbit using 



= 4 1 



2E n Jl 

Gy?M 2 ' 



where M = Mi + M 2 and fj, = MiM 2 /M, 

GM l M 2 
an - 2^T' 

r n = (1 + e„)a n , 
and 



T n = 



47r 2 a n 
GM 



(5) 



(6) 



(T) 



We can then evolve the binary system according to our 
assumptions. The results of applying this procedure to the 
encounter between a 9 Mq star and a 1 Mq star at a perias- 
tron of 10 Rq are shown in Figure |7| The effect of the tidal 
interaction between the stars is to circularise their orbit, so 
the eccentricity, semimajor axis and orbital period decrease 
rapidly until, after w 115 yr, the apastron is approximately 
equal to the periastron and the peraistron itself begins to de- 
crease. At this point we assume that the system will rapidly 
merge. 

We determined the binary grind-down times for all of 
our encounters and found that they ranged from a few hun- 
dred hours for encounters with periastrons of 4 Rq to « 100 
yr for encounters between 9 Mq stars and 1 Mq stars with 
periastrons of 10 Rq or encounters between two 9 Mq stars 
wit h periastro n s of 15 Rq . 

IMardlind (Il995al) and iMardlind il995bT) and also 
iKocha nck (1992) studied the evolution of tidal capture bi- 
naries using numerical schemes in which the cores of the 
stars were represented by point masses and the polytropic 
envelope of the larger star was decomposed into spherical 
harmonics. The point masses and envelope components were 
allowed to exchange energy under the assumption that the 
overall interaction was adiabatic. This allowed the authors 
to integrate systems similar to those formed in our calcula- 
tions for many thousands of orbits. Mardling found that, for 
very eccentric orbits (as the orbits of our captured systems 
are, since they were captured from parabolae), the evolution 
of the binary separation was unpredictable but that, in gen- 
eral, very e ccent ric binaries did not grind down to mergers. 
iKochanekl (Il992l) found that binaries formed by encounters 
at a few stellar radii were often scattered into wide orbits on 
subsequent periastron passages and that on, average, such 
binaries spent most of their time as soft binaries and were 
therefore vulnerable to disruption. 

There is a third, more likely outcome of tidal capture. 
If the binary grind-down times t gr ind and final separations 
a/ estimated by our method above are realistic, the tidal 
power Ludai deposited in the envelopes may be estimated 
from 



GM1M2 



(8) 



(4) 



2^ grind 

The tidal luminosity of the encounter between a 9 Mq 
star and a 1 Mq star is then ~ 10 Lq, far in excess of 
the luminosity of a main-sequence 9 Mq star. It is there- 
fore likely that the tidal interaction will cause the star(s) 
to expand and enter a phase of common-envelope evolution. 
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Figure 7. Evolution of the orbital eccentricity, period and semimajor axis in the captured system formed by the encounter of a 9 Mq 
star and a 1 Mq star with a periastron of 10 Rq. 



3.3 Common envelope evolution 

If our captured systems do grind down tidally, the likely re- 
sult is a common-envelope system in which the stellar cores 
inspiral inside a common envelope derived from the stars' 
original envelopes. The fate of a common envelope system 
hinges on whether the energy released during the inspiral is 
sufficient to unbind the envelope before the cores merge. 

If the common-envelope phase is assumed to expel the 
envelope and if the energy released derives purely from the 
initial orbital energy of the inspiralling cores, one may de- 
rive an expression for the final separation of the cores, qf as 
a func tion of their initial separation di (e.g. iDewi fc Taurisl 
<2C 



.... f 2(M C + M env )Menv , Mi{M c + M env ] 

at = Ivliivic — 1 

\ aXR at 

where Mi is the mass of the intruder star, M c is the core 
mass of the larger star, M env is the mass of the larger star's 
envelope, a is a parameter expressing the efficiency with 
which the orbital energy of the cores is used to unbind the 
envelope and A is a geometrical factor accounting for the en- 
velope's structure. Formally, a may take any positive value; 
q = 1 implies that all the orbital binding energy is required 
to eject the common envelope, and a > 1 implies that some 
additional source of energy is required. To determine the 
outcome of common-envelope evolution for a given system, 
the product qA can be treated as a parameter. 

If the separation between the cores becomes so small 
that one of them fills its Roche lobe (note that, since these 
systems are not corotating, the use of Roche geometry is 
an approximation, but it proves to be adequate), the cores 
merge and energy deposition in the envelope ceases. For two 
cores i and j of mass rm and m,j at a s eparation a, th e size 
of the Roche lobe of core i is given by feggletonl llQSah 

0A9(m,i/mj)ia 



r' R = , V " t ' / " t ^ " -. (10) 

0.6(mi/mj) 3 + ln(l + (rrii/mj)a) 

We continuously estimated the Roche lobe size for both 
cores from their masses and separation and using Equation 
1101 as the simulation evolved. We compared the estimated 
Roche lobe size of each core to that core's physical size. The 
cores of the initial model stars were defined by the hydrogen- 



burning regions in the models' centres and all SPH particles 
within this region were tagged at the start of each simula- 
tion. We then estimated the sizes of the stellar cores from the 
positions of all the core particles, assuming that each core 
neither lost nor gained particles and that the cores were 
spherical. At the point in time where the radius of either 
core calculated in this way first exceeded the Roche lobe 
size of that core, an animation of the simulation was exam- 
ined to confirm that mass transfer between the cores was 
in progress. In both cases, we used the encounters with a 
periastron of 2 Rq . In encounters between 9 Mq stars and 1 
Mq stars, core Roche-lobe overflow occurs at a separation 
of ~ 1.2 Rq and in encounters between two 9 Mq stars, it 
occurs at a separation of « 1.7 Rq. Taking the initial separa- 
tions for both types of encounter to be the relevant capture 
(9) radii (as would be the case following the tidal grind-down 
of a captured system) , we may now use Equation [5] to de- 
termine the minimum value of aX that would result in each 
system ejecting its envelope before undergoing core merger. 

For encounters between 9 Mq stars and 1 Mq stars, we 
took the core, envelope and intruder masses to be 2 Mq, 7 
Mq and 1 Mq respectively. We found that we would require 
extreme values of aX (> 20) to eject the envelope before 
core merger occurs in the case of an inspiral from the cap- 
ture radius for these encounters, 11 Rq. We repeated the 
analysis for encounters between two 9 Mq stars, taking the 
core, envelope and intruder masses to be 2 Mq, 14 Mq and 
2 Mq respectively. In the case of inspiral from a radius of 18 
Rq (the capture radius for these encounters) , ejection of the 
envelope requires even more extreme values of aX, > 100. 
The value of A for the initial 9 Mq star models was ~ 7. This 
value can be expected to increase as the 9 Mq star expands, 
but we would probably still require values of a in excess 
of unity to expel the envelope, implying that some source 
of energy other than the orbital energy of the inspiralling 
cores would be required. Such a scenario is very difficult to 
study in the energy-conserving model of common-envelope 
evolution, but another prescriptio n for studying common- 
env elope systems was suggest ed by Neleman^^^Toud (2004) 
iNelemans fc Toud (l2004h and lNelemans fc Tout] ho03) 
developed an approach to common envelope evolution based 
on explicit conservation of angular momentum. They assume 
that the fractional quantity of angular momentum carried 
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away by ejection of the envelope is proportional to the frac- 
tional mass loss, so that 

AJ = AM _ M env 

J ~ 7 M ~ J M C + M env ' [ ' 

One may then derive an expression for a/ as a function of 
a; analogous to Equation [3] 



(M env + M c )Mj 
M c Mi 



M c + Mi 



M env + M c + M t 



1-7 



M e „„ + M c + Mi 



(12) 



In their study of white dwar f binaries as the p r oducts 
of common envelope evolution, iNelemans fc Tout! J2005T) 
concluded that the origins of all the systems they exam- 
ined could be explained using the above expressions if 
1.5 < 7 < 1.75. We applied their analysis to our stellar 
encounters, again assuming that the grind-down of a cap- 
tured binary resulted in common-envelope evolution where 
the initial separation between the two cores corresponded 
to the appropriate capture radius. We found that the 
analysis of Nelemans and Tout (2004) predicts that, under 
these conditions, encounters between 9 Mq stars and 1 Mq 
stars will result in expulsion of the envelope for 7 > 1.65 
and that encounters between two 9 Mq stars will result 
in envelope expulsion for any value of 7 in Nelemans and 
Tout's preferred range. In both cases, the final separations 
of the systems were very sensitive to the chosen value of 7. 

The results of the calculations performed in this 
section are inconclusive. We have applied two different 
common-envelope formalisms to encounters between 9 Mq 
stars and 1 Mq stars and between two 9 Mq stars and 
found that the traditional explicitly energy-conserving 
formalism implies that envelope ejection is unlikely and 
that these encounters probably result in mergers, whereas 
th e explicitly angu l ar-mo mentum-conserving formalism 
of INelemans fc Touti (120051) implies that envelope ejection 
is possible in both cases and that the stellar cores will 
not merge. This essentially means that the orbits of these 
common-envelope systems have sufficient angular mo- 
mentum, but insufficient energy, to expel their envelopes. 
Settling the question of whether or not the systems merge 
rests on whether some additional source of energy can 
be tapped, so that the envelope can be e xpelled while 
conserving energy and angular momentum. ISokerl <l2004h 
suggests that angular momentum deposition causes the 
envelope to swell and triggering the formation of dust. 
The resulting dust-driven wind taps the stellar luminosity, 
which provides the 'extra' source of energy required to eject 
the envelope. In the absence of an additional energy source, 
the result of common-envelope evolution would presumably 
be that a small quantity of material would be expelled from 
the systems, carrying all the orbital angular momentum, 
although it is not obvious how this might come about. 

We have only considered massive stars of a single mass. 
On the main sequence, the fraction of the stellar mass 
made up by the core is an increasing function of stellar 
mass. We repeated our common envelope analysis for an 
encounter between a 50 Mq star, with a 35 Mq core and 
a 15 M Q envelope llKippenhahn fc Weigerdll994T) . and a 1 
M© intruder. We found that a merger could be avoided and 
the stellar envelope ejected in common-envelope evolution 



following binary grind-down for qA > 5 and for any value 
of 7 in the range preferred by Nelemans and Tout. This 
result is counterintuitive at first sight, but follows from 
the fact that the core of the 50 Mq star has such a large 
mass, so that there is a large quantity of orbital energy 
available with which to expel the envelope, and also that 
the envelope has a relatively small mass and binding energy. 
Sim ilar conclusions hav e been reached by ot her authors, 
e.g. iTerman et"afl (Il994l) . iTerman etUl i 19951) . who found 
that a 24 Mq supergiant star entering a common-envelope 
phase with a 1.4 Mq neutron star would lose its envelope. 



3.4 Encounters where rp < (Ri + R2) 

In this section, we discuss encounters in which the perias- 
tron distance on the first approach is less than the sum of 
the stellar radii, so that the stars collide immediately. The 
dissipative interaction between the two stars ejects mass 
from the system and the cores of the stars spiral inwards. 

We integrated the encounters between 9 Mq stars and 
1 Mq stars with periastrons of 1, 2 and 3 Rq and those 
between two 9 Mq stars with periastrons of 1, 2, 3 and 
4 Rq until the cores of the stars merged (doing this for 
runs with larger periastrons was not possible, since the 
stars were captured onto relatively wide orbits which would 
require prohibitive quantities of CPU time to integrate). 
In all the simulations which we ran to core merger, the 
mergers were very rapid (~ 15 — 80 hours in the encounters 
between 9 Mq stars and 1 Mq stars and « 12 — 150 hours 
in the encounters between two 9 Mq stars). The total 
quantities of mass ejected by the mergers were small in 
comparison to the sum of the progenitor masses (~ 0.1 
Mq — 0.3 Mq in encounters between 9 Mq stars and 1 Mq 
stars and « 0.4 Mq — 0.9 Mq in the encounters between 
two 9 Mq stars). These results are largely consistent 
with our common-envelope calculations in the sense that 
the calculations presented in Section 3.3 suggest that 
encounters between 9 Mq stars and 1 Mq stars with these 



periastrons should not expel the envelope of the 9 Mq 
star for any v alue of a A or for a ny va lue of 7 considered 
reasonable bv INelemans fc Touti (l2004h . Envelope ejection 
in encounters between two 9 Mq stars with periastrons 
^ 4 Rq is forbidden by the energy-conserving analysis, 
although the 7-formalism suggests that envelope ejection 
may occur for values of 7 > 1.6. Note, however, that we 
are not here suggesting that common-envelope evolution 
is generally amenable to study by SPH simulations, since 
integration for many dynamical times would be required 
in the majority of cases. Apart from the prohibitively long 
CPU times required, a more accurate artificial viscosity 
prescription would probably be required, since significant 
momentum transport by viscosity may occur over a large 
number of orbits. In our SPH calculations, mergers occur 
in only a few tens of dynamical times and almost all 
momentum transfer or transport is by gravitational torques 
or the ram-pressure interaction between the stars. It would 
also be difficult to include other possible sources of energy, 
such as dust-driven winds. Our SPH calculations implicitly 
assume that no such sources of energy are available and 
hence that a does not exceed unity. 

To determine the probable subsequent evolution of the 
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Figure 8. Histogram showing to which star material in 0.5 Mq 
shells within the collision product of the encounter between a 9 
Mq star and a 1 Mq star with a periastron of 2 Rq originally 
belonged. Hatched material originated in the 9 Mq star and white 
material in the 1 Mq star. 

products of these mergers, we examined the distribution 
of the material from the progenitor stars in the collision 
products. In Figure |H] we divide the collision product of 
the encounter between a 9 Mq star and a 1 Mq star with 
a periastron of 2 Rq into spherical 0.5 Mq shells and show 
what fraction of the material in each shell originated from 
the 1 Mq star and what fraction came from the 9 Mq star. 

The 1 Mq star has survived the merger largely intact 
and occupies the centre of the merger product. Since the 1 
Mq star is considerably denser and cooler than the 9 Mq 
star, it is not surprising that it has sunk in this way. In 
the corresponding encounter between two 9 Mq stars, we 
observed that the collision product was well mixed after 
~ 30 h ours. We conclude, in c omm on with [ Lo mbardi et alJ 
l|2002l) , ISandauist et alJ |l99Tt> and ISills et all (|2002l ) fand 
contrary to iBenz fc Hills! <ll987l) ) that little hydrodynamic 
mixing occurs during collisions of main sequence stars of 
unequal mass, although it is possible that mixing will occur 
when the merger product is allowed to approach thermal 
equilibrium, since the core will then become convective. 

In Figure [5] we plot the fraction by mass of helium 
against enclosed mass in the collision product and com- 
pare the results from our SPH calculation from results 
derived using t h e entr opy-sorting technique suggested by 
lLombardi et alJ (120021) . The two methods evidently agree 
closely, both suggesting that the merger remnant will have 
a core with a helium-rich centre surrounded by a layer rich 
in fresh hydrogen, in turn surrounded by a layer enriched 
in helium. If the merger remnant eventually comes to 
approximate a normal main sequence 10 Mq star, it will 
have a core mass of ~ 3 Mq. The combined core masses of 
the progenitor stars is » 2.3 Mq, so the hydrogen-burning 
region will need to grow before the remnant settles onto the 
main sequence. The innermost 3 Mq of the merger remnant 
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Figure 9. Plot of the helium fraction by mass against en- 
closed mass within the collision product of the encounter be- 
tween a 9 Mq star and a 1 Mq star with a periastron of 2 
Rq. We compare the results of the SPH calculation with those 
derived using the entropy-sorting method of 'L ombardi et alJ 
( 2002 ) . The dashed plot was generated by inserting our initial 
stellar models into the Make Me A Star package available at 
http : //faculty . vassar . edu/lombardi/mmas/. 

has a mean helium fraction of ~ 70%, whereas the core of 
the original 9 Mq star had a core helium fraction of ~ 90%. 
The duration of the main-sequence phase for the model 
9 Mq star generated using the ATON code was 3.0 x 10 7 
yr. We used the model at an age of 2.8 x 10 7 yr in our 
SPH calculations. If the merger product eventually settles 
back onto the main sequence, it should continue to burn 
hydrogen for w 6 x 10 6 yr, whereas the original 9 Mq star 
had enough hydrogen left to burn for only another w 2 x 10 6 
yr. The merger may then extend the hydrogen-burning 
lifetime of the 9 M star by w 4 x 10 6 yr, w 13% of its 
total main-sequence lifetime. This is an extreme case of 
rejuvenation. If the 9 Mq star were younger and thus had 
a less helium-enriched core, the effect of the merger on the 
core helium fraction would clearly be smaller. The merger 
of two evolved stars, as in our mergers of two 9 Mq stars, 
also clearly does not result in a remnant with a much 
extended hydrogen-burning phase, since the remnant core 
is formed from two already highly helium-enriched cores. 



3.5 Massive stars with rapidly— rotating cores - 
hypernova progenitors? 

Hypernovae (HNe) are a class of unusually luminous 
supernovae, speci fically Type Ic even ts with luminosities 
>fewxl0 51 erg ijNomoto et alJ Il999f) . Recently, a con- 
nection between HNe and long-duration Gamma Ray 
Bursts (GRBs) (all references to GRBs in this paper 
relate to long-perio d burs ts) has been posited (see, e.g. 
iPodsiadlowski et alJ (120041) . suggesting that they may be 
different facets of the same type of event, which is itself a 
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high-energy subset of Type Ic supernovae. HNe/GRB mod- 
els involving stellar progenitors require extreme rotation 
rates to ensure that the collapsing stellar core forms a disk 
around its central black hole. Although young massive stars 
are often rapidly rotating, there are numerous means by 
which they can lose their angular momentum during their 
main-sequence lifetimes (e.g. via winds). 

Mergers are an ideal candidate for producing rapidly- 
rotating massive stars, since all the orbital angular 
momentum of the system must reside either in the merger 
product or in the small amount of mass ejected in the 
collision. Mergers in which little mass is lost therefore result 
in remnants with very large angular momenta. 

For the core of a massive star to form a HN/GRB, the 
specific angular momentum of material at the edge of the 
core before collapse begins must exceed the specific angular 
momentum of material on the last stable orbit around the 
black hole, so that the whole core cannot collapse into the 
hole, giving j > y/QGM cor Jc » 2 x 10 16 (M core /2M Q ) erg 
s g^ JPodsiadlowski et al]l2004f) . The merger of two 9 M 
stars may therefore result in a hypernova progenitor if the 
core of the product is rotating sufficiently rapidly. 

We examined the remnant of the encounter between 
two 9 Mq stars with a peri astron of 2 Rq to see if it fulfils 
iPodsiadlowski et alJ l2004l) 's criterion. Defining the core 
radius as the radius containing 2 M , the core radius of the 
remnant is ~ 0.8 Rq and the specific angular momentum at 
this radius is ~ fifty times that required to prevent the core 
of the collision product collapsing directly into a central 
black hole. We were only able to analyse the collision 
remnant « 100 hr after impact and angular momentum 
transport is still underway within the remnant, largely 
by means of spiral density inhomogeneities. However, the 
remnan t can lose 98% of its ang ular momentum and still 
satisfy iPodsiadlowski et al.ll20o'3) 's criterion. 



4 DISCUSSION 

4.1 Comparison with runs using point mass 
impactors 

As well as the high-resolution calculations described in Sec- 
tion 3, we conducted two simulations of encounters between 
9 Mq stars and 1 Mq stars with a periastron of 2 R in 
which we represented the 1 Mq star by a point mass with 
gravitational smoothing lengths of 0.5 and 0.05 Rq, chosen 
so that the point mass smoothing length would be com- 
parable to the SPH smoothing lengths in the 9 Mq star's 
envelope and core respectively. We compare these runs with 
the pure-SPH run with a periastron of 2 Rq. 

The main difference between the pure-SPH and point- 
mass runs was that there was greater dissipation in the pure- 
SPH runs. The SPH intruder star experiences ram pressure, 
which the point mass does not, and thus loses kinetic en- 
ergy more rapidly during its passag e through the 9 Mq star. 
This phenomenon is discussed bv lBenz fc Hilla al992f) . who 
also performed collision simulations using point-mass and 
SPH impactors, and can be seen clearly in Figure [T01 which 
contrasts the time-evolution of the core separation in the 
pure-SPH and point-mass runs, and in Figure [Til where we 
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Figure 10. Comparison of the evolution of the core separations 
with time in the pure-SPH and point mass encounters between 9 
Mq stars and 1 Mq stars with periastrons of 2 Rq. Dot— dashed 
line: pure-SPH calculation, thick solid line: point mass calculation 
with gravitational smoothing length of 0.05 Rq, thin solid line: 
point mass calculation with gravitational smoothing length of 0.5 
Rq. 

compare the quantities of rotational kinetic energy trans- 
ferred to the 9 Mq star on the first few periastron passages 
in each run. The smaller dissipation in the point-mass sim- 
ulations results in merger timescales which are longer than 
those in the pure-SPH calculations, but are still only of the 
order of tens of hours. 

The mass-loss evolution of these calculations is shown 
m Figure [HI The mass-loss in the point-mass simulations 
is greater because gas particles in these runs are strongly 
accelerated in the deep potential wells of the point masses. 
A smaller gravitational smoothing length produces a deeper 
well and thus higher a cceleration (behaviou r of this kind was 
discussed at length bv lDavies et alJ (11995^ 1. When compar- 
ing results of runs using SPH and point-mass intruders, the 
point masses will always eject more mass in this way than an 
equal-mass SPH intruder, unless their gravitational smooth- 
ing lengths are considerably larger than the SPH intruder 
itself. 

We conclude from these results that replacing an SPH 
impactor with a point mass in stellar collision calculations 
produces qualitatively similar results, but is unsuitable if 
one wants to study the energetics of the encounters quanti- 
tatively. Point mass impactors are more suita ble for studies 
of enc o unters involving compa ct objects, e.g. iTerman et all 
Jl995l) . lLavcock fc Sillsl J2005f) . 



4.2 Merger products as hypernova progenitors 

We showed in Section 3 that the cores of collision remnant 
from encounters between two 9 Mq stars are rotat- 
ing sufficiently rapidly to fulfil the criterion outlined 
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Figure 11. Evolution of the rotational kinetic energy of the 9 Mq 
star during the first few periastron passes in the pure-SPH and 
point mass encounters between 9 Mq stars and 1 Mq stars with a 
periastron of 2 Rq. Dot— dashed line: pure-SPH calculation, thick 
solid line: point mass calculation with gravitational smoothing 
length of 0.05 Rq , thin solid line: point mass calculation with 
gravitational smoothing length of 0.5 Rq. 
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Figure 12. Comparison of the evolution of the mass loss in the 
pure-SPH and point mass encounters between 9 Mq stars and 1 
Mq stars with a periastron of 2 Rq. Dot— dashed line: pure-SPH 
calculation, thick solid line: point mass calculation with gravita- 
tional smoothing length of 0.05 Rq, thin solid line: point mass 
calculation with gravitational smoothing length of 0.5 Rq. 



by IPodsiadlowski et alJ i2004f) . making them plausible 
candidates for HNe/GRB progenitors. We now consider 
whether the cluster environments in which massive stars 
are typically found permit collision rates sufficiently large 
to account for a significant fraction of HNe. 

The stellar densities in globular cluster cores often 
reach 10 4 pc -3 , although the massive stars in such systems 
have long since exploded. If we consider younger systems 
which still possess most of their high-mass stars, we find 
that, for example, the Or ion Nebula Cluster has a stella r 
density of ~ 10 4 pc -3 iMcCaughrean fc Stauffeil I1994T) . 
the putative Galactic starb urst cluster NGC 36 03 has a 
mass density of ~ 10 5 pc -3 llHofmann et al.l|l995t) and the 
Arches cluster near t he Galactic Centre h as a mass density 
of > 5 x 10 5 pc" 3 l)Figer fc Kiml l2002h . If we take 10 4 
pc -3 as being representative of typical clusters containing 
massive stars, we can use Equation 1141 to estimate the 
fraction of massive main-sequence stars in these systems 
that will undergo collisions. Taking these stars to have radii 
of 10 Rq, masses of 10 Mq and taking the cluster velocity 
dispersion to be 10 km s , we find that the encounter 
timescale is ~ 10 10 yr. If the massive stars in question have 
main-sequence lifetimes of ~ 10 7 yr and they spend their 
entire lives in a cluster environment, ~ 0.1% of them will 
suffer collisions in clusters with number densities of 10 4 
pc -3 . 

The t otal rate of Type Ib/c supe rnovae is ~ 10 -3 yr -1 
galaxy -1 (IPodsiadlowski et alJ feOO^ I. Podsiadlowski et al 
estimate that 5% of observed SNe Ib/c are HNe. Taking 
into account the fact that HNe are brighter than average 
SNe Ib/c, they estimate that the HNe rate is ~ 10 -5 yr -1 
galaxy -1 . The total rate of all core-collapse supernovae is 
~ 7 x 10 -3 yr -1 galaxy -1 , implying that ~ 0.1% of all 
massive stars explode as HNe. We therefore find that the 
fraction of massive stars likely to suffer collisions in stellar 
systems with number densities typical of young massive 
clusters is of the same order as the fraction of massive stars 
that explode as hypernovae. 

To explode as a Type Ic supernova, the merger remnant 
must first lose its envelope. Winds are an obvious means 
of achieving this and common-envelope evolution may also 
expel a significant fraction of the envelope of the merger 
product during the merger. 



4.3 Ultimate fate of collision remnants 

We have so far discussed only single collisions between main 
sequence stars, such as occur infrequently in clusters of high 
stellar density. We showed in Section 4.2 that in clusters with 
densities of ~ 10 4 pc -3 , ~ 0.1% of massive stars are likely to 
suffer collisions. This fraction increases linearly with stellar 
density, reaching 100% at a density of ~ 10 7 pc -3 . The ex- 
treme stellar densities (~ 10 s pc -3 ) though t to exist in some 
young stellar systems jBonnell et al]IT998l) bring with them 
the possibility of multiple collisions. To understand the fate 
of an encounter involving the product of an earlier merger, 
the evolution of merger remnants must be considered. 

The gravitationally-focused cross section of stars of 
mass M and radius R in a stellar cluster with velocity dis- 
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persion Voc is given by 

If Voo is small (a value of 10 km s _1 is assumed throughout 
this paper) , the second term inside the brackets of Equation 
| 13l d ominates and an encounter timescale can be derived 
jFabian et al"lll975l) 

---""(^rarGs)"* 

( V °° A yr, (14) 
VlOkms" 1 / ' V ; 

where n is the stellar number density of the cluster. In this 
paper we have neglected the possibility that the presence 
of primordial binaries or the formation of tidal-capture 
binaries will probably decrease the effective encounter 
timescales, since binaries have m uch larger encounter 
cross- sections than single stars. iPortegies Zwart et alJ 
( 1999) suggest that binary-formation may increase collision 
rates by a factor of ~ 10. 

In the encounters between 9 Mq stars and 1 Mq stars, 
almost all of the kinetic energy from the 1 Mq star was 
deposited in the merger product's envelope. The quantity 
of energy deposited in the envelope of the 9 Mq star in the 
encounter with a periastron of 2 Rq was ~ 10 7 Lq yr. The 
luminosity of a 10 Mq star is 10 3 — 10 4 Lq, so this energy 
will be radiated away on a timescale of 10 3 - 10 4 yr. The 
thermal timescale of the remnant is ~ 10 4 yr. The merger 
remnant is therefore likely to swell in a manner similar 
to a star leaving the main sequence and evolving into a 
giant. If the luminosity of the merger product remains 
const ant, the post-main sequ e nce tr ack for a 10 Mq star 
(e.g. iKippenhahn fc Weigerd lll994f) ') suggests that the 
remnant will swell by a factor of ~ 100 to a radius of ~ 500 
Rq. After ~ 10 4 yr, the collision remnant will have radiated 
away the excess energy deposited in its envelope and will 
shrink back to the size of a main-sequence 10 Mq star. 
The collision cross section of the bloated remnant, (Tbioat 
is ~ 100 times that of the remnant on the main sequence, 
oms and the time for which the remnant remains bloated, 
tbioat is a few x 10 4 years, several hundred to a thousand 
times shorter than its likely survival time if it contracts onto 
the main sequence undisturbed, tias- Let the probability 
of the remnant being struck a second time while in its 
bloated state be Pbioat oc (Tbloattbloat an d the probability 
of the remnant being struck again after contracting to the 
main sequence be Pms oc (TMstMS- It follows from the 
rough cancellation of the factors by which auoat > cms 
and by which tbioat < tias that, for any given remnant, 
Pbioat ~ Pms- In stellar clusters where the probability of 
repeated collisions is low, this relation is unimportant since 
Pbioat and Pms are both very small, but it becomes more 
interesting if stars are likely to suffer more than on collision. 

We repeated the common-envelope analyses from 
Section 3.3 to determine the likely outcome of an encounter 
between a bloated collision remnant and a main-sequence 
star, again taking the remnant's capture radius to be 1000 
Rq. For a A > 7 or for any value of 7 within Nelemans and 
Tout's preferred range, we found that such an encounter 
was unlikely to lead to a merger, but instead would eject 
the remnant's envelope, leaving a binary composed of the 



cores of the merger remnant and of the second intruder 
star. In stellar clusters where t enc ~ tbioat, Pbioat becomes 
significant and it follows that roughly half of collision 
remnants will be struck a second time in their bloated state 
and disrupted, while roughly half will be struck again after 
contracting onto the main sequence. The results of Section 
3.5 suggest that those remnants surviving to complete their 
main-sequence phase may explode as hypernovae. 

If we take the remnant's capture radius to be 1000 
Rq, its mass to be 10 Mq and set the encounter timescale 
in Equation 1141 to the duration of the bloated state, we 
find that the number density at which this occurs is a 
fewxlO 7 pc -3 (inserting instead the main-sequence radius 
and lifetime gives a similar result). 

In denser stellar systems (n ~ 10 s pc -3 ) the population 
of remnants contracting onto the main sequence is depleted, 
since the majority of remnants suffer subsequent collisions 
while bloated and never reach the main sequence. Since 
the encounter timescale is a few times shorter than the 
thermal timescale of the product of the first collision, 
merger remnants will be struck several times while swollen. 
If the timescale on which the cores inspiral into the centre 
of the object and merger is comparable to or longer than 
the collision time, at any one instant, the envelope of the 
remnant may contain three or four cores, forming a bag of 
cores. This state of affairs would lead to a complex few- 
body encounter, possibly expelling the lower-mass core(s) 
while common-envelope evolution ejects the envelope, 
leaving a binary composed of the two highest-mass cores, 
or alternatively lead ing to the merger of two of the cores 
jFreeeau et alJl2004l) . If, instead, the timescale for the cores 
to inspiral and merge is shorter than the collision timescale, 
each new impactor will probably merge with the object 
undisturbed by interactions with other impactors. The 
evolution of such systems will be the topic of a subsequent 
papery 

IPorteeies Zwart fc McMi llan| (|200 2[) an d more r ecently 
iFreitae et all (l2005bK and iFreitag et alJ J2005d) have 
suggested that clusters of still higher number densities 
(10 9 — 10 11 pc -3 , reached during core collapse) may host a 
runaway merger process in which tens of merger events lead 
to the formation of an object of mass ~ 10 3 Mq which sub- 
sequently explodes and leaves behind an intermediate-mass 
black hole. We do not think this scenario is very likely in 
practice. At these stellar densities, the encounter timescales 
are much shorter than the thermal timescale of the remnant 
of the initial merger, so collision products are pummelled 
by tens of impacts before they get the opportunity to swell 
on a thermal timescale. However, they will still swell on a 
dynamical timsecale, so later mergers may eject significant 
quantities of mass. In addition, massive stars have very 
powerful winds which drive mass-loss rates that can reach 
10~ 4 Mq yr" 1 even for 'ordinary' OB or Wolf-Rayet stars 
Ide Jager et all 1988). The wind from a star with a mass of 
~ 1000 Mq is likely to be stronger still. It is therefore likely 
that the runaway buildup of mass becomes self-limiting, 
so that the rate at which the merger product gains mass 
by mergers is comparable to the rate at which mass is 
being ejected. Even if a very massive object does form, 
once collisions cease and the remnant approaches thermal 
equilibrium, it will swell to a vast radius and develop a 
very pronounced core/envelope structure, so that most 
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Case Number density Number of collisions per object and fate of collision products 
(pc-3) 



(0) 


n < 10 4 


Collisions rare - stars evolve normally. 


(1) 


10 < n < 10 


0.1 — lO/o of stars will sutler single collisions. Products go on to explode as rapidly— rotating 
supernovae (RSN) or hypernovae if sufficiently massive (> 20 Mq). 


(2) 


n ~ 10 7 


Most stars suffer 1 — 2 collisions. About half of products likely to be struck again while 

swollen and to lose their envelopes, leaving a hard binary. The other half suffer a second 

collision on the main sequence and go on to explode as rapidly-rotating supernovae (RSN) or hypernovae 

if sufficiently massive (> 20 Mq). 


(3) 


n ~ 10 8 


All stars likely to suffer several collisions. Collision remnants likely to be struck two or more 
times while in their swollen state. Interaction of the cores with each other and with the envelope 
likely to eject the lowest-mass cores as well as the envelope, leaving highest-mass cores in a hard 
binary. 


(4) 


10 9 < n < 10 11 


All stars likely to be involved in collisions. Collision timescale much shorter than remnants' thermal 
timcscales so remnants have no opportunity to swell. Multiple mergers may build up an object 
of mass ~ 100 Mq after which mass buildup probably self-limited by winds. 



Table 2. Likely numbers of collisions and fates of collision products in stellar environments of increasing number density. Note that 
any given star cluster may pass through some or all of these regimes as it evolves and may exhibit several of these regimes at any given 
moment in time. 



of its mass is concentrated in a very small volume at its 
centre, the rest being distributed in a tenuous envelope. 
This structure, together with the tremendous luminosity 
that the core of such an object must have, will almost 
certainly lead to insta bilities driving mass loss r ates of at 
least ~ 10" 3 M yr _1 jGlatzel fc Kiriakididll993l) . 

We summarise the probable fates of collision remnants 
in stellar environments of increasing number density in Ta- 
ble|3] In Figure we explain schematically the sequence of 
stellar encounters corresponding to each row in Tabled The 
stellar number density in clusters, particularly in the cores 
of dense clusters, can increase dramatically on timescales 
comparable to or shorter than massive stars' main-sequence 
lifetimes due to relaxation-driven core-collapse and this 
evolution can be influenced by st ellar collisions, sinc e 
collisions can accelerate core collapse (iFreitag et al.ll2005al ). 
Conversely, the expulsion of gas be feedback from young 
stars can decrease the core stellar number density (e.g. 
iKrouoa et alJ l|200lll )In addition, the number density in a 
given cluster clearly varies spatially, generally being higher 
in the core than at, for example, the half-mass radius. 
Table and Figure EH are intended to show the expected 
fate of stars in environments with the number densities 
indicated. A given cluster may pass through some or all of 
these density regimes as it evolves and may exhibit some or 
all of these regimes at any given moment in time. 



5 CONCLUSIONS 

We have studied encounters between 9 Mq stars and 1 Mq 
stars and between two 9 M stars in stellar systems with 
velocity dispersions of ~ 10 km s~ . The capture radius 
for encounters between 9 Mq stars and 1 M stars with 
this velocity dispersion is w 11 R and that for encounters 
between two 9 Mq stars is « 18 R . 

Using high- and very high-resolution SPH simulations, 
we showed that encounters with small periastrons (r p 3 
Rq in encounters between 9 Mq stars and 1 Mq stars, 
and r p ^ 4 Rq in encounters between two 9 Mq stars) 



resulted in mergers taking place on timescales of tens to 
hundreds of hours. By comparing results obtained at two 
very different numerical resolutions, we ascertained that 
such short merger times are not an artefact of numerical 
resolution. 

Encounters with larger periastrons led to captured 
systems whose orbits would probably shrink over timescales 
of hundreds of hours to hundreds of years, leading to 
a phase of common envelope evolution. We studied the 
evolution of the common envelope systems resulting 
from the grind-down of the captured systems using the 
standard energy-conserving common-envelope formal- 
ism and the angular-m omentum-conserving method of 
iNelemans fc Touq 12004). The energy-conserving formalism 
suggested that common-envelope encounters resulting 
from the grind-down of a captured system would result 
in a merger unless extreme values of the product aX were 
invoked - > 20 in encounters between 9 Mq stars and 1 
M stars, and > 100 in encounters between two 9 M 
stars. The results for the same encounters derived using the 
angular-momentum-conserving formalism were different 
and very sensitive to the chosen value of the 7 parameter. 
We again considered common-envelope evolution following 
the grind-down of a captured system. In encounters be- 
tween 9 Mq stars and 1 Mq stars, we found that values 
°f 7 si 1-65 resulted in envelope-ejection without core 
merger. In encounters between two 9 Mq stars, any value 
within Nelemans and Tout's preferred range resulted in 
envelope-ejection without core merger. We also found 
that more massive stars (~ 50 Mq), with their relatively 
more massive cores and less massive envelopes, are more 
vulnerable to envelope expulsion. 

We integrated our SPH simulations of the encounters 
between 9 Mq stars and 1 Mq stars with periastrons ^ 3 
Rq and those between two 9 Mq stars with periastrons 
^ 4 R until the stellar cores merged. These mergers 
ejected small quantities of mass, and, in mergers of 9 Mq 
stars and 1 Mq stars, the smaller star spiralled largely 
intact into the core of the larger star. The 9 Mq models 
used in our calculations had sufficient core hydrogen for 
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Figure 13. Likely fate of the products of collisions in environments with the number densities given in Table l2l 'RSN' indicates that 
these collision remnants will explode as rapidly-rotating supernovae (and those with masses of > 20 Mq possibly as hypernovae). As 
explained in the text, any stellar cluster may pass through some or all of these density regimes as it evolves and may exhibit some or all 
regimes at any given moment in time. 



k 2 x 10 yr of main sequence evolution. Merging these 
objects with 1 Mq stars delivered enough fresh hydrogen 
to the cores of the merger remnants to allow them to burn 
for another ~ 6 x 10 6 yr. We conclude that merger with a 
lower-mass object may extend the duration of a massive 
star's main-sequence phase by up to ~ 10%, although 
such an extension is likely to be smaller in reality, since we 
considered massive stars close to terminal age. 

The products of our collisions were rapidly rotating. 
The merger of two 9 Mq stars with a periastron of 2 R© 
produced an object with a core rotating sufficiently rapidly 
to make it a candidate hypernova/gam ma ray burst progen- 
itor ac cording to the criterion given in IPodsiadlowski et alJ 
(2004). Common-envelope evolution may expel some or 
all of the remnants' envelopes, so that their subsequent 
supernova/hypernova explosions are of Type Ib/c. Taking 



the typical stellar clusters in which massive stars are found 
to have stellar number densities of ~ 10 4 pc -3 , we found 
that ~ 0.1% of massive stars will suffer collisions, although 
only remnants with masses > 20 Mq are massive enough 
to eve ntually explode as hypernovae. IPodsiadlowski et alJ 
(2004) estimate that ~ 0.1% of massive stars explode as 
hypernovae. We showed that merger remnants possess 
sufficient core specific angular momentum after their for- 
mation to make them possible HNe/GRB progenitors, that 
the typical encounter timescales in massive stellar systems 
imply a collision rate similar to the hypernova rate and 
suggest that common-envelope evolution and winds may 
expel the envelopes of collision remnants. We conclude that 
the products of massive stellar collisions may contribute 
significantly to the observed hypernova rate. 

We studied further interactions of our collision prod- 
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ucts with main-sequence stars. The collision remnants swell 
on a thermal timescale, reaching radii large enough that 
an encounter with another star would result in a phase of 
common envelope evolution which expels the remnant's 
envelope. However, once the excess energy has been radiated 
away, an undisturbed remnant contracts (also on a thermal 
timescale) and becomes a main-sequence star. In stellar 
clusters with number densities of ~ 10 7 pc -3 , about half of 
the collision products suffer a second collision back on the 
main sequence while the other half are struck whilst swollen 
and lose their envelopes, leaving a hard binary. In systems 
with number densities of ~ 10 s pc -3 , the collision product 
will suffer multiple further collisions whilst swollen, leading 
to a system in which several (i.e.>two) cores orbit within a 
common envelope. Interactions between the cores may eject 
all but the two most massive and drive off the envelope, 
leaving a hard binary, or alternatively lead to a merger of 
two of the cores. We will investigate this problem in detail 
in a subsequent paper. In denser systems (n > 10 9 pc -3 ) 
the remnant of the first collis io n will be struc k many times. 
IPortegies Zwart & McMillan! <l2002h and iFreitag et ail 
(2005aj) and IFreitag et al.l ll2005b!) have suggested that 
such a runaway merger process will result in the formation 
of a very massive (~ 1000 Mq) object which eventually 
explodes and forms an intermediate-mass black hole. We do 
not think this scenario is likely as hydrodynamic swelling 
of the remnant, making it vulnerable to mass-ejection 
in subsequent collisions, and strong stellar winds would 
probably make the mass buildup self limiting when the 
merger product acquires a mass of a few hundred solar 
masses. In addition, the collision re mnant may b ecome 
'transparent' to subsequent imapctors .Even 
if a ~ 1000 Mq object does form, we suggest that winds 
and instabilities are likely to result in mass losses > 10 3 Mq 
yr _ . By the time the object reaches the end of its main 
sequence phase, it is likely to have lost most of its mass 
and its supernova explosion will not produce an unusually 
massive black hole. 
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